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Abstract 

The Anderson-Hubbard Hamiltonian at half-filling is investigated within dynamical mean-field theory at zero 
temperature. The local density of states is calculated by taking the geometric and arithmetic mean, respectively. 
The non-magnetic ground state phase diagrams obtained within the different averaging schemes are compared. 
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The metal insulator transitions (MIT) due to elec- 
tronic interactions (Mott-Hubbard MIT) [1] and due 
to impurity scattering (Anderson localization) [2] are 
subtle quantum mechanical phenomena which require 
non-perturbative investigation schemes. In this respect 
the dynamical mean-field theory (DMFT) [3] is very 
useful. However, it cannot describe the physics of An- 
derson localization if the local disorder is included by 
taking the arithmetic average over the disorder [4] . Re- 
cently it was proposed to use the geometric average 
over the disorder to include also the Anderson MIT 
[5,6,7]. In this contribution we compare the entire non- 
magnetic ground state phase diagrams of correlated, 
disordered electrons at half filled lattice as obtained by 
taking the geometric and the arithmetic averages over 
the disorder within the DMFT [8]. 

We study the system described by a single-orbital 
Anderson-Hubbard model 

Hah = -t 22 a l<r a 3" + ^2 t%nia + n *T n a> (!) 
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where t > is the amplitude for hopping between near- 
est neighbors, U is the on-site repulsion, n; CT = a\ a a ia 
is the local electron number operator, (aJ CT ) is the 
annihilation (creation) operator of an electron with 
spin a, and the local ionic energies e< are independent 
random variables. We assume a continuous probabil- 
ity distribution for e;, i.e., V(ei) = 0(A/2 — [ej[)/A, 
with B as the step function. Here, A is a measure of 
the disorder strength. 

This model is solved within DMFT by mapping it [3] 
onto an ensemble of effective single-impurity Anderson 
Hamiltonians with different e;: 

#siam = 5^(ei — n)a\ a aia + Umtmx (2) 

ker kcr 

Here \i — U/2 is the chemical potential corresponding 
to a half-filled band, and Vk and ek are the hybridiza- 
tion matrix element and the dispersion relation of the 
auxiliary bath fermions Cko- , respectively. 

For each ionic energy we calculate the local Green 
function G(u), from which we can obtain either 
the geometrically averaged local density of states 
(LDOS) pgeom(w) = exp [(ln/9j(w))] or the arith- 
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Fig. 1. Non-magnetic ground state phase diagram of the An- 
derson-Hubbard model at half-filling as calculated by DMFT 
with the geometrically averaged local density of states. I — 
correlated, disordered metal, II — coexistence regime, III — 
crossover regime, IV - Mott insulator, V — Anderson insula- 
tor. At the dashed line all states within the Hubbard subbands 
become localized. W is a bare energy band-width. 

metically averaged LDOS p a rith(w) = {pi(ui)), where 
Pi(uj) = —lmG(Lu,ei)/iv, and (Oi) = J deiP(ei)0(ei) 
is the arithmetic mean of Oi. The lattice Green func- 
tion is given by the corresponding Hilbert transform 
as G(u>) = J du>' p a (uj') / (lo — lo'), where the subscript 
a stands for either "geom" or "arith". The local 
self-energy £(u>) is determined from the k- integrated 
Dyson equation = u — t](uj) — 1/G(u) where 

the hybridization function r;(w) is defined as rj(uo) = 
Ek l^kl 2 / (uj - e k ). The self-consistent DMFT equa- 
tions are closed through the Hilbert transform G(lu) = 
f deJVo(e)/ [u — e — ~£(lu)], which relates the local 
Green function for a given lattice to the self-energy; 
here iVu(e) is the non-interacting DOS. 

The Anderson-Hubbard model (1) is solved for a 
semi-elliptic DOS, No(e) = 4VT - 4e 2 /tt. Then t)(lu) = 
G(w)/16. The DMFT equations are solved at zero tem- 
perature by the numerical renormalization group tech- 
nique [9]. For numerical integrations we use discrete 
values of e; selected according to the Gauss-Legendre 
algorithm. The number of e; levels depends on A and 
is adjusted to obtain smooth spectral functions [4]. 

The main results are presented in Fig. 1 and 2, where 
the non-magnetic ground state phase diagrams of the 
Anderson-Hubbard model (1) are plotted. At weak dis- 
order both averaging schemes are seen to lead to very 
similar results for the Mott-Hubbard transition from 
the paramagnetic correlated, disordered metal (I) to 
the paramagnetic, disordered Mott insulator (IV). For 
A < 2 we find hysteresis and a coexistence regime 
(II). Above A ~ 2 there is a crossover regime (III). 
At strong disorder the DMFT with geometrically aver- 
aged LDOS describes an Anderson transition between 
a correlated, disordered metallic and insulating phase, 
respectively (V) [10]. Since neither the paramagnetic 
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Fig. 2. Non— magnetic ground state phase diagram of the An- 
derson—Hubbard model at half-filling as calculated by DMFT 
with the arithmetically averaged local density of states. Labels 
I — IV are the same as in Fig. 1. 

Mott nor the Anderson insulator is characterized by a 
broken symmetry they are continously connected, i.e., 
the phases IV and V in Fig. 1 are connected by a con- 
tinuous line which does not cross a metallic phase. Ob- 
viously, the Anderson transition is missed within the 
DMFT supplied by the arithmetic averaging. 

In summary, although the geometric averaging pro- 
cedure does not capture all aspects of Anderson lo- 
calization [6,11] it provides valuable new insights into 
MITs in correlated and disordered electron systems, 
which are not obtained by taking the arithmetic disor- 
der average [4]. 
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